Comparative plastome analysis and taxonomic classification of snow lotus species (Saussurea, Asteraceae) in Central Asia and Southern Siberia.

Four species of Saussurea, namely S. involucrata, S. orgaadayi, S. bogedaensis, and S. dorogostaiskii, are known as the "snow lotus," which are used as traditional medicines in China (Xinjiang), Kyrgyzstan, Mongolia, and Russia (Southern Siberia). These species are threatened globally, because of illegal harvesting and climate change. Furthermore, the taxonomic classification and identification of these threatened species remain unclear owing to limited research. The misidentification of medicinal species can sometimes be harmful to health. Therefore, the phylogenetic and genomic features of these species need to be confirmed. In this study, we sequenced five complete chloroplast genomes and seven nuclear ITS regions of four snow lotus species and other Saussurea species. We further explored their genetic variety, selective pressure at the sequence level, and phylogenetic relationships using the chloroplast genome, nuclear partial DNA sequences, and morphological features. Plastome of the snow lotus species has a conserved structure and gene content similar to most Saussurea species. Two intergenic regions (ndhJ-ndhK and ndhD-psaC) show significantly high diversity among chloroplast regions. Thus, ITS and these markers are suitable for identifying snow lotus species. In addition, we characterized 43 simple sequence repeats that may be useful in future population genetic studies. Analysis of the selection signatures identified three genes (rpoA, ndhB, and ycf2) that underwent positive selection. These genes may play important roles in the adaptation of the snow lotus species to alpine environments. S. dorogostaiskii is close to S. baicalensis and exhibits slightly different adaptation from others. The taxonomic position of the snow lotus species, confirmed by morphological and molecular evidence, is as follows: (i) S. involucrata has been excluded from the Mongolian flora due to misidentification as S. orgaadayi or S. bogedaensis for a long time; (ii) S. dorogostaiskii belongs to section Pycnocephala subgenus Saussurea, whereas other the snow lotus species belong to section Amphilaena subgenus Amphilaena; and (iii) S. krasnoborovii is synonymous of S. dorogostaiskii. This study clarified the speciation and lineage diversification of the snow lotus species in Central Asia and Southern Siberia.

Complete Chloroplast Genomes of Saussurea katochaete, Saussurea superba, and Saussurea stella: Genome Structures and Comparative and Phylogenetic Analyses.

Saussurea plants are widely distributed in Asia and Europe; however, their complex phylogenetic relationships have led to many difficulties in phylogenetic studies and interspecific identification. In this study, we assembled, annotated, and analyzed the chloroplast genomes of three Saussurea plants: Saussurea katochaete, Saussurea superba, and Saussurea stella. The results showed that the full-length sequences of the three Saussurea plants were 152,561 bp, 151,452 bp, and 152,293 bp, respectively, which represent the typical quadripartite structure, and the genomes were relatively conserved. The gene annotation results showed that the chloroplast genomes of S. katochaete, S. superba, and S. stella were annotated with 128, 124, and 127 unique genes, respectively, which included 83, 80, and 83 protein-coding genes (PCGs), respectively, 37, 36, and 36 tRNA genes, respectively, and 8 rRNA genes. Moreover, 46, 45, and 43 SSR loci, respectively, and nine highly variable regions (rpl32-trnL-UAG, rpl32, ndhF-rpl32, ycf1, trnC-GCA-petN, trnC-GCA, rpcL, psbE-petL, and rpl16-trnG-UUG) were identified and could be used as potential molecular markers for population identification and phylogenetic study of Saussurea plants. Phylogenetic analyses strongly support the sisterhood of S. katochaete with S. superba and S. stella, and are all clustered with S. depsagensis, S. inversa, S. medusa, and S. gossipihora, of which S. gossipiphora is most closely related. Additionally, the phylogenetic results indicate a high frequency of differentiation among different species of Saussurea plants, and many different species or genera are morphologically very different from each other, which may be related to certain genetic material in the chloroplasts. This study provides an important reference for the identification of Saussurea plants and studies their evolution and phylogenetics.

Metabolomics-based exploration the response mechanisms of Saussurea involucrata leaves under different levels of low temperature stress.

BACKGROUND: Saussurea involucrata (Sik.) is alpine plant that have developed special adaptive mechanisms to resist adverse environmental conditions such as low temperature chilling during long-term adaptation and evolution. Exploring the changes of its metabolites under different temperature stresses is helpful to gain insight into its cold stress tolerance. METHODS: Ultra-performance liquid chromatography and tandem mass spectrometry were used to analyze the metabolites in the leaves of Sik. under low different temperature stress conditions. RESULTS: A total of 753 metabolites were identified, and 360 different metabolites were identified according to the Kyoto Encyclopedia of Genes and Genomes (KEGG) involved in the biosynthesis of secondary metabolites and amino acids and sugars. Sucrose and trehalose synthesis, glycolysis, tricarboxylic acid cycle, pentose phosphate pathway, glutamic acid-mediated proline biosynthesis, purine metabolism, amino acid metabolism, phenylpropane synthesis pathway metabolites all respond to low temperature stress. Under cold stress conditions, carbohydrates in Sik. leaves accumulate first than under freezing conditions, and the lower the temperature under freezing conditions, the less amino acids accumulate, while the phenolic substances increase. The expression of various substances in LPE and LPC increased more than 10-fold after low temperature stress compared with the control, but the content of LPE and LPC substances decreased after cold adaptation. In addition, purines and phenolics decreased and amino acids accumulated significantly under freezing conditions. CONCLUSION: The metabolic network of Sik. leaves under different low temperature stress conditions was proposed, which provided a reference for further exploration of the metabolic mechanism related to low temperature stress tolerance of Sik.

Integrated analysis of smRNAome, transcriptome, and degradome data to decipher microRNAs regulating costunolide biosynthesis in Saussurea lappa.

Saussurea lappa (S. lappa) has been known to synthesize medicinally important, costunolide. Due to its immense therapeutic importance, understanding of regulatory mechanism associated with its biosynthesis is crucial. The identification of genes and transcription factors (TFs) in S. lappa, created a clear picture of costunolide biosynthesis pathways. Further to understand the regulation of costunolide biosynthesis by miRNAs, an integrated study of transcriptome, miRNAs, and degradome was performed. Identified candidate miRNAs and associated feed-forward loops (FFLs) illustrates their regulatory role in secondary metabolite biosynthesis. Small RNA and degradome sequencing were performed for leaf and root tissues to determine miRNAs-targets pairs. A total of 711 and 525 such targets were obtained for novel and known miRNAs respectively. This data was used to generate costunolide-specific miRNA-TF-gene interactome to perform systematic analyses through graph theoretical approach. Interestingly, miR171c.1 and sla-miR121 were identified as key regulators to connect and co-regulate both mevalonate and sesquiterpenoid pathways to bio-synthesize costunolide. Tissue-specific FFLs were identified to be involved in costunolide biosynthesis which further suggests the evolutionary co-relation of root-specific networks in synthesis of secondary metabolites in addition to leaf-specific networks. This integrative approach allowed us to determine candidate miRNAs and associated tissue-specific motifs involved in the diversification of secondary metabolites. MiRNAs identified in present study can provide alternatives for bioengineering tool to enhance the synthesis of costunolide and other secondary metabolites in S. lappa.

RD29A promoter constitutively drives a rice Hd3a expression to promote early-flowering in Saussurea involucrate Kar. et Kir. ex Maxim.

S. involucratae, an endemic and endangered plant, is a valuable and traditional Chinese medicinal herb. In order to control the flowering time of S. involucratae, we used the well-known stress inducible RD29A promoter to drive Hd3a (a FT ortholog from rice) expression in S. involucratae. Unexpectedly, the majority of regenerated buds in RD29A::Hd3a transgenic lines (S-RH) produced flowers in tissue culture stage under normal growth (25 ± 2 °C) condition. Their flowering time was not further influenced by salt treatment. Hd3a in S-RH was strongly expressed in MS media supplemented with or without 50 mM NaCl. RD29A::GUS transgenic experiments further revealed that RD29A constitutively promoted GUS expression in both S. involucrate and halophyte Thellungiella halophile, in contrast to glycophic plants Oryza sativa L. 'Zhonghua 11', in which its expression was up-regulated by cold, salinity, and drought stress. The results supported the hypothesis that RD29A promoter activity is inducible in stress-sensitive plants, but constitutive in stress-tolerant ones. Importantly, S-RH plants produced pollen grains and seeds under normal conditions. Additionally, we found that OsLEA3-1::Hd3a and HSP18.2::Hd3a could not promote S. involucrate to flower under either normal conditions or abiotic stresses. Taken together, we demonstrated the potential of RD29A::Hd3a might be served as a feasible approach in breeding S. involucrate under normal condition.

Saussurea involucrata PIP2;4 improves growth and drought tolerance in Nicotiana tabacum by increasing stomatal density and sensitivity.

Aquaporins, the major facilitators of water transport across membranes, are involved in growth and development and adaptation to drought stress in plants. In this study, a plasma membrane intrinsic protein (SiPIP2;4) was cloned from Saussurea involucrata, a cold-tolerant hardy herb. The expression of SiPIP2;4 increased the stomatal density and sensitivity of tobacco (Nicotiana tabacum), thus, affecting the plant's growth and resistance to the diverse water environment. The higher stomatal density under well-watered conditions effectively promoted the photosynthetic rate, which led to the rapid growth of transgenic lines. The stomata in the transgenic lines responded more sensitively to the vapor pressure deficit than the wild-type under different levels of ambient humidity. Their stomatal apertures positively correlated with the ambient humidity. Under drought conditions, the overexpression of SiPIP2;4 promoted rapid stomatal closure, reduced water dissipation, and enhanced drought tolerance. These results indicate that SiPIP2;4 regulates the density and sensitivity of plant stomata, thus, playing an important role in balancing plant growth and stress tolerance. This suggests that SiPIP2;4 has the potential to serve as a genetic resource for crop improvement.

Comparative plastomes and phylogenetic analysis of seven Korean endemic Saussurea (Asteraceae).

BACKGROUND: Saussurea is one of the most species-rich genera in the Cardueae, Asteraceae. There are approximately 40 Saussurea species distributed in Korea, with nearly 40% of them endemics. Infrageneric relationships remain uncertain due to insufficient resolutions and low statistical support. In this study, we sequenced the plastid genomes of five Korean endemic Saussurea (S. albifolia, S. calcicola, S. diamantica, S. grandicapitula, and S. seoulensis), and comparative analyses including two other endemics (S. chabyoungsanica and S. polylepis) were conducted. RESULTS: The plastomes of Korean endemics were highly conserved in gene content, order, and numbers. Exceptionally, S. diamantica had mitochondrial DNA sequences including two tRNAs in SSC region. There were no significant differences of the type and numbers of SSRs among the seven Korean endemics except in S. seoulensis. Nine mutation hotspots with high nucleotide diversity value (Pi > 0.0033) were identified, and phylogenetic analysis suggested that those Korean endemic species most likely evolved several times from diverse lineages within the genus. Moreover, molecular dating estimated that the Korean endemic species diverged since the late Miocene. CONCLUSIONS: This study provides insight into understanding the plastome evolution and evolutionary relationships of highly complex species of Saussurea in Korean peninsula.

Selection and Validation of Reference Genes for RT-qPCR Analysis in Tibetan Medicinal Plant Saussurea Laniceps Callus under Abiotic Stresses and Hormone Treatments.

Real-time quantitative PCR (RT-qPCR) is an important technique for studying gene expression analysis, but accurate and reliable results depend on the use of a stable reference gene. This study proposes to test the expression stability of candidate reference genes in the callus of Saussurea laniceps, a unique Tibetan medicinal plant. Based on the S. laniceps callus transcriptome, eleven candidate reference genes, including TUA2, TUB3, TUB8, TIF3B1, TIF3H1, ELF5A, PP2AA2, UEV1D, UBL5, UBC36, and SKIP1), were validated for RT-qPCR normalization in the callus under abiotic stress (salt, cold, and UV) and hormone treatments (abscisic acid, MeJA, and salicylic acid). The stability of the candidate genes was evaluated in all the samples of S. laniceps. Comprehensive analysis of all samples showed that the best reference genes were UBC36 and UBL5. ELF5A and TIF3B1 were ranked as the most stable genes in the sample sets under abiotic stress. For hormone stimulation, UBC36 and TIF3H1 genes had the best stability. This study provides useful guidelines and a starting point for reference gene selection for expression analysis using RT-qPCR techniques in S. laniceps.

Macroevolutionary pattern of Saussurea (Asteraceae) provides insights into the drivers of radiating diversification.

Evolutionary radiations have intrigued biologists for more than a century, yet our understanding of the drivers of radiating diversification is still limited. We investigate the roles of environmental and species-intrinsic factors in driving the rapid radiation of Saussurea (Asteraceae) by deploying a number of palaeoenvironment-, diversity- and trait-dependent models, as well as ecological distribution data. We show that three main clades of Saussurea began to diversify in the Miocene almost simultaneously, with increasing diversification rates (DRs) negatively dependent on palaeotemperature but not dependent on species diversity. Our trait-dependent models detect some adaptive morphological innovations associated with DR shifts, while indicating additional unobserved traits are also likely driving diversification. Accounting for ecological niche data, we further reveal that accelerations in DRs are correlated with niche breadth and the size of species' range. Our results point out a macroevolutionary scenario where both adaptive morphological evolution and ecological opportunities provided by palaeoenvironmental fluctuations triggered an exceptionally radiating diversification. Our study highlights the importance of integrating phylogenomic, morphological, ecological and model-based approaches to illustrate evolutionary dynamics of lineages in biodiversity hotspots.

Saussurea involucrata (Snow Lotus) ICE1 and ICE2 Orthologues Involved in Regulating Cold Stress Tolerance in Transgenic Arabidopsis.

As with other environmental stresses, cold stress limits plant growth, geographical distribution, and agricultural productivity. CBF/DREB (CRT-binding factors/DRE-binding proteins) regulate tolerance to cold/freezing stress across plant species. ICE (inducer of CBF expression) is regarded as the upstream inducer of CBF expression and plays a crucial role as a main regulator of cold acclimation. Snow lotus (Saussurea involucrata) is a well-known traditional Chinese herb. This herb is known to have greater tolerance to cold/freezing stress compared to other plants. According to transcriptome datasets, two putative ICE homologous genes, SiICE1 and SiICE2, were identified in snow lotus. The predicted SiICE1 cDNA contains an ORF of 1506 bp, encoding a protein of 501 amino acids, whereas SiICE2 cDNA has an ORF of 1482 bp, coding for a protein of 493 amino acids. Sequence alignment and structure analysis show SiICE1 and SiICE2 possess a S-rich motif at the N-terminal region, while the conserved ZIP-bHLH domain and ACT domain are at the C-terminus. Both SiICE1 and SiICE2 transcripts were cold-inducible. Subcellular localization and yeast one-hybrid assays revealed that SiICE1 and SiICE2 are transcriptional regulators. Overexpression of SiICE1 (35S::SiICE1) and SiICE2 (35S::SiICE2) in transgenic Arabidopsis increased the cold tolerance. In addition, the expression patterns of downstream stress-related genes, CBF1, CBF2, CBF3, COR15A, COR47, and KIN1, were up-regulated when compared to the wild type. These results thus provide evidence that SiICE1 and SiICE2 function in cold acclimation and this cold/freezing tolerance may be regulated through a CBF-controlling pathway.

Genetic diversity and structure of Saussurea polylepis (Asteraceae) on continental islands of Korea: Implications for conservation strategies and management.

Saussurea polylepis Nakai is an herbaceous perennial endemic to Korea and is highly restricted to several continental islands in the southwestern part of the Korean Peninsula. Given its very narrow geographical distribution, it is more vulnerable to anthropogenic activities and global climate changes than more widely distributed species. Despite the need for comprehensive genetic information for conservation and management, no such population genetic studies of S. polylepis have been conducted. In this study, genetic diversity and population structure were evaluated for 97 individuals from 5 populations (Gwanmaedo, Gageodo, Hongdo, Heusando, and Uido) using 19 polymorphic microsatellites. The populations were separated by a distance of 20-90 km. We found moderate levels of genetic diversity in S. polylepis (Ho = 0.42, He = 0.43). This may be due to long lifespans, outcrossing, and gene flow, despite its narrow range. High levels of gene flow (Nm = 1.76, mean Fst = 0.09), especially from wind-dispersed seeds, would contribute to low levels of genetic differentiation among populations. However, the small population size and reduced number of individuals in the reproductive phase of S. polylepis can be a major threat leading to inbreeding depression and genetic diversity loss. Bayesian cluster analysis revealed three significant structures at K = 3, consistent with DAPC and UPGMA. It is thought that sea level rise after the last glacial maximum may have acted as a geographical barrier, limiting the gene flow that would lead to distinct population structures. We proposed the Heuksando population, which is the largest island inhabited by S. polylepis, as a source population because of its large population size and high genetic diversity. Four management units (Gwanmaedo, Gageodo, Hongdo-Heuksando, and Uido) were suggested for conservation considering population size, genetic diversity, population structure, unique alleles, and geographical location (e.g., proximity).

Identification and characterization of SlbHLH, SlDof and SlWRKY transcription factors interacting with SlDPD gene involved in costunolide biosynthesis in Saussurea lappa.

The genes involved in costunolide biosynthesis in Saussurea lappa have been identified recently by our lab. However, the study of transcriptional regulators of these genes was lacking for better opportunities for engineering the pharmacologically important biosynthetic pathway. Therefore, we cloned the promoter region of diphosphomevalonate decarboxylase gene (DPD) and analyzed its cis-acting regulatory elements to reveal the potential transcription factor (TF) binding sites for Dof, bHLH and WRKY family proteins in the gene promoter. The transcriptome study approach followed by the hidden Markov model based search, digital gene expression, co-expression network analysis, conserved domain properties and evolutionary analyses were carried out to screen out seven putative TFs for the DPD-TF interaction studies. Yeast one-hybrid assays were performed and three TFs were reported, namely, SlDOF2, SlbHLH3 and SlWRKY2 from Dof, bHLH and WRKY families, respectively that interacted positively with the DPD gene of the costunolide biosynthetic pathway. The tissue specific relative gene expression studies also supported the linked co-expression of the gene and its interacting TFs The present report will improve the understanding of transcriptional regulation pattern of costunolide biosynthetic pathway.

Comparative transcriptomics reveals candidate transcription factors involved in costunolide biosynthesis in medicinal plant-Saussurea lappa.

Costunolides, an important sesquiterpene lactone (STL) isolated from Saussurea lappa, are the major pharmaceutical ingredient of various drug formulations. Identification of the genes and transcriptional regulation of costunolide biosynthesis pathway in S. lappa will propose alternatives for engineering enhanced metabolite biosynthesis in plant. Here, we aimed to unravel the transcription factors (TFs) regulating the costunolide biosynthesis. Comparative transcriptome analysis of root and leaf tissues and transcripts were annotated using various in silico tools. Putative transcription factors were identified using PlantTFDB and TF- gene co-expression network was generated followed by clustering using module based analysis to observe their coordinated behaviour. The module 1 was found to be significant based on its enrichment with major pathway genes. Further, promoter cloning determined the cis acting elements in costunolide synthase (SlCOS1) gene which catalyses the final key step of costunolide biosynthesis. Bioinformatics tools were employed to predict the cis regulatory elements, leading to the identification of MYB family of TFs as an interacting partner of SlCOS1 gene. The present study is the pioneer attempt for TF prediction and elucidation of their regulatory role in costunolide synthesis. This will help in future metabolic engineering of the pharmaceutically important STLs and their yield improvement.

Overexpression of Saussurea involucrata dehydrin gene SiDHN promotes cold and drought tolerance in transgenic tomato plants.

Dehydrins are late embryogenesis abundant proteins that help regulate abiotic stress responses in plants. Overexpression of the Saussurea involucrata dehydrin gene SiDHN has previously been shown to improve water-use efficiency and enhance cold and drought tolerance of transgenic tobacco. To understand the mechanism by which SiDHN exerts its protective function, we transformed the SiDHN gene into tomato plants (Solanum lycopersicum L.) and assessed their response to abiotic stress. We observed that in response to stresses, the SiDHN transgenic tomato plants had increased contents of chlorophyll a and b, carotenoid and relative water content compared with wild-type plants. They also had higher maximal photochemical efficiency of photosystem II and accumulated more proline and soluble sugar. Compared to those wild-type plants, malondialdehyde content and relative electron leakage in transgenic plants were not significantly increased, and H2O2 and O2- contents in transgenic tomato plants were significantly decreased. We further observed that the production of stress-related antioxidant enzymes, including superoxide dismutase, ascorbate peroxidase, peroxidase, and catalase, as well as pyrroline-5-carboxylate synthetase and lipid transfer protein 1, were up-regulated in the transgenic plants under cold and drought stress. Based on these observations, we conclude that overexpression of SiDHN gene can promote cold and drought tolerance of transgenic tomato plants by inhibiting cell membrane damage, protecting chloroplasts, and enhancing the reactive oxygen species scavenging capacity. The finding can be beneficial for the application of SiDHN gene in improving crop tolerance to abiotic stress and oxidative damage.

A novel cold-regulated protein isolated from Saussurea involucrata confers cold and drought tolerance in transgenic tobacco (Nicotiana tabacum).

Adverse environmental conditions, such as cold and drought, can inhibit plant growth, development, and productivity. The isolation and characterization of stress response genes from stress-tolerant plants can provide a better understanding of the underlying adaptive mechanisms. In this study, a novel cold-regulated gene, SikCOR413PM1, was isolated from Saussurea involucrata Kar. et Kir., which is a plant that survives at the high altitudes and in the low temperatures of alpine slopes in northwestern China. SikCOR413PM1 was induced in response to cold and drought in S. involucrata, and phylogenetic analysis revealed that the gene groups with a COR gene encoding a COR413PM protein family member. Subcellular localization of a SikCOR413PM1-green fluorescent fusion protein showed that SikCOR413PM1 was localized to the plasma membrane. A transgenic tobacco (Nicotiana tabacum) system was employed to investigate the possible role of SikCOR413PM1 in cold and drought tolerance. Analyses of growth, germination and survival rates, relative water content, malondialdehyde content, relative electrolyte leakage, and maximal photochemical efficiency of photosystem II showed that transgenic tobacco plants expressing SikCOR413PM1 were more tolerant to cold and drought stresses than WT plants. SikCOR413PM1 overexpression was also accompanied by constitutive activation of NtDREB1 and NtDREB3, two cold-responsive transcription factor genes, and NtERD10A and NtERD10B, two cold-induced genes. The expression levels of downstream transcription factor genes NtDREB3, NtERD10C, NtERD10D, and NtLEA5 were also induced in SikCOR413PM1-expressing transgenic plants under drought conditions. Our results suggest that the overexpression of SikCOR413PM1 induces changes in tobacco plants, and facilitates enhanced tolerance to cold and drought stresses.

Phylogeny, origin and dispersal of Saussurea (Asteraceae) based on chloroplast genome data.

Saussurea is one of the largest genera of the tribe Cardueae of Asteraceae, comprising about 460 species from the Northern Hemisphere with most species distributed in QTPss and adjacent areas. Here, we established a well-supported phylogenetic framework for Saussurea based on whole chloroplast genomes of 136 taxa plus 16 additional taxa of Cardueae using Bayesian inference and Maximum Likelihood. Our phylogenetic results are inconsistent with previous subgeneric classifications of Saussurea. We nearly completely delimited subgen. Eriocoryne, and found that subgen. Theodorea, subgen. Saussurea section Laguranthera and Rosulascentes are closely related to each other. Based on our phylogenetic results, we performed biogeographic analyses and inferred that the genus Saussurea arose during early-middle Miocene within the Hengduan Mountains. We expect that landscape heterogeneity within the QTPss and adjacent areas, such as the Hengduan Mountains, played an important role in the evolution of Saussurea. Following its evolutionary origin, the genus underwent rapid diversification in situs and dispersed northwards in several migrational patterns. Both continuous uplift of the QTPss and adjacent areas as well as global cooling since mid-Miocene probably led to geographic expansion and diffusion of Saussurea, with the latter, in particular, resulting in the northward dispersal.

Plastome phylogenomics of Saussurea (Asteraceae: Cardueae).

BACKGROUND: Saussurea DC. is one of the largest and most morphologically heterogeneous genera in Asteraceae. The relationships within Saussurea have been poorly resolved, probably due an early, rapid radiation. To examine plastome evolution and resolve backbone relationships within Saussurea, we sequenced the complete plastomes of 17 species representing all four subgenera. RESULTS: All Saussurea plastomes shared the gene content and structure of most Asteraceae plastomes. Molecular evolutionary analysis showed most of the plastid protein-coding genes have been under purifying selection. Phylogenomic analyses of 20 Saussurea plastomes that alternatively included nucleotide or amino acid sequences of all protein-coding genes, vs. the nucleotide sequence of the entire plastome, supported the monophyly of Saussurea and identified three clades within it. Three of the four traditional subgenera were recovered as paraphyletic. Seven plastome regions were identified as containing the highest nucleotide variability. CONCLUSIONS: Our analyses reveal both the structural conservatism and power of the plastome for resolving relationships in congeneric taxa. It is very likely that differences in topology among data sets is due primarily to differences in numbers of parsimony-informative characters. Our study demonstrates that the current taxonomy of Saussurea is likely based at least partly on convergent morphological character states. Greater taxon sampling will be necessary to explore character evolution and biogeography in the genus. Our results here provide helpful insight into which loci will provide the most phylogenetic signal in Saussurea and Cardueae.

Elucidating genes involved in sesquiterpenoid and flavonoid biosynthetic pathways in Saussurea lappa by de novo leaf transcriptome analysis.

Saussurea lappa (family Asteraceae) possesses immense pharmacological potential mainly due to the presence of sesquiterpene lactones. In spite of its medicinal importance, S. lappa has been poorly explored at the molecular level. We initiated leaf transcriptome sequencing of S. lappa using the illumina highseq 2000 platform and generated 62,039,614 raw reads. Trinity assembler generated 122,434 contigs with an N50 value of 1053 bp. The assembled transcripts were compared against the non-redundant protein database at NCBI. The Blast2GO analysis assigned gene ontology (GO) terms, categorized into molecular functions (3132), biological processes (4477) and cellular components (1.927). Using KEGG, around 476 contigs were assigned to 39 pathways. For secondary metabolic pathways, we identified transcripts encoding genes involved in sesquiterpenoid and flavonoid biosynthesis. Relatively low number of transcripts were also found encoding for genes involved in the alkaloid pathway. Our data will contribute to functional genomics and metabolic engineering studies in this plant.

New Saussurea (Asteraceae) species from Bogeda Mountain, eastern Tianshan, China, and inference of its evolutionary history and medical usage.

In this study, Saussurea bogedaensis Yu-J. Wang & Jie Chen, a new species from Bogeda Mountain in the eastern part of the Tianshan Mountains, is described and discussed based on evidence in terms of both morphological and genetic data. S. bogedaensis is morphologically similar to S. involucrata, which is distributed in the western part of the Tianshan Mountains, and it is well known because of its beauty, rarity, and medicinal value. The new species is also similar to S. orgaadayi, which is distributed in the nearby Altai Mountains. Our genetic data support the close relationships among these three species. According to their allopathic distributions, we suggest that these three species are derived from the same ancestor but that they differentiated after reaching their current range. In addition, we propose that the new species might serve as an alternative to S. involucrata in medicine due to their very high similarity. However, this species appears to be rare because we only found six mature individuals in the field despite extensive investigations.

De Novo Transcriptome Sequencing and the Hypothetical Cold Response Mode of Saussurea involucrata in Extreme Cold Environments.

Saussurea involucrata grows in high mountain areas covered by snow throughout the year. The temperature of this habitat can change drastically in one day. To gain a better understanding of the cold response signaling pathways and molecular metabolic reactions involved in cold stress tolerance, genome-wide transcriptional analyses were performed using RNA-Seq technologies. A total of 199,758 transcripts were assembled, producing 138,540 unigenes with 46.8 Gb clean data. Overall, 184,416 (92.32%) transcripts were successfully annotated. The 365 transcription factors identified (292 unigenes) belonged to 49 transcription factor families associated with cold stress responses. A total of 343 transcripts on the signal transduction (132 upregulated and 212 downregulated in at least any one of the conditions) were strongly affected by cold temperature, such as the CBL-interacting serine/threonine-protein kinase (CIPKs), receptor-like protein kinases, and protein kinases. The circadian rhythm pathway was activated by cold adaptation, which was necessary to endure the severe temperature changes within a day. There were 346 differentially expressed genes (DEGs) related to transport, of which 138 were upregulated and 22 were downregulated in at least any one of the conditions. Under cold stress conditions, transcriptional regulation, molecular transport, and signal transduction were involved in the adaptation to low temperature in S. involucrata. These findings contribute to our understanding of the adaptation of plants to harsh environments and the survival traits of S. involucrata. In addition, the present study provides insight into the molecular mechanisms of chilling and freezing tolerance.